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Abstract. 



Recurrence coefficients of semi-classical orthogonal polynomials (or- 
thogonal polynomials related to a weight function w such that w' /w 
is a rational function) are shown to be solutions of non linear dif- 
ferential equations with respect to a well-chosen parameter, accord- 
ing to principles established by D.& G. Chudnovsky. Examples are 
given. For instance, the recurrence coefficients in an+iPn+i{x) = 
xpn{x) — anPn-i{x) of the orthogonal polynomials related to the 
weight exp(— — tx"^) on M satisfy 4a^ari = (3a^ -|- 2ta^ — n)(a^ + 
2ta^ + n) , and satisfies a Painleve Piv equation. 



1. Introduction: measures and recurrence coefficients of orthogonal 
polynomials. 

Let {pn}o° the set of orthonormal polynomials related to some measure d/i on its 
support 5": 



Pnix)Pmix) dp{x) = 6m,n- (1) 

S 

The most remarkable property of the p^s is the recurrence relation joining them: 

an+lPn+l{x) = (X - bn)Pnix) - a„p„_i(x). (2) 

An often encountered problem in applied and numerical mathematics as well as in physics 
is then to relate the coefficients and bn of (2) to properties of the measure d/j,. 



For instance, interesting solid-state Hamiltonian operators submitted to the so- 
called "recursion method" (or Lanczos method) show a tri-diagonal matrix representation. 
Investigation of spectral properties of the operator is then equivalent to investigating the 
measure of (1) from the recurrence coefficients of (2) [GaCL] [Hay] [HayN] [LaG] [LiMu] 
[OW]. 
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The study of special partition functions in statistical physics and quantum physics 
leads to relations which can be translated as properties of particular recurrence coefficients. 
Much important work is currently done on this subject [Bes] [BIZ] [Fokl] [Fok2] [Fral] 
[Fra2] [GrMl] [GrM2] [HH] [KvM] [LW] [Mo] [Y] [Zu]. 

Numerical implementation of spectral methods and quadrature formulas needs 
accurate determination of recurrence coefficients for various measures. This appears in 
the survey [Gau] and in some recent papers as [Chin] and [CIS] (see also the references in 
[BeR]). 

Quite a number of theoretical studies have appeared on this problem of relating 
properties of the orthogonality measure to the recurrence coefficients, especially to their 
asymptotic behaviour. See at least the books [Chi] , [Frl] , [VA] and the survey papers [Lub] 
and [GFOPCF]. 

To give just a taste of the matter, the asymptotic behaviour of the recurrence 
coefficients associated to diJ,{x) — Ix]'' exp{—x'^)dx on x gM appears in an amazing number 
of fields: 

1. This extension of Hermite polynomials is studied by Shohat [Sho], using methods of 
Laguerre [Lag]. Later on, Freud [Fr2] rediscovered Shohat 's formulas (see (31) with 
t = 0) and proved that ~ (n/12)^/'^ when n oo. Much more has been done 
since then [Lub] [Mag2] [Mag3] [Nev] [GFOPCF] [Nev2], it has been shown that a 
behaviour dii{x) ~ exp(— for x — > ±oo implies an and hn ~ constant n^/" for 
large n. 

2. Similar exponential weights were encountered in solid-state and statistical physics, 
where the same asymptotic connections have been used (sometimes after clever 
guesswork) [LiMu] [OW]. 

3. These extensions of Hermite and Laguerre polynomials also appear in numerical 
quadrature methods intended to solve Boltzmann and Fokker-Planck equations [CIS] , 
where they are called "speed", "bimodc" and "Druyvestcyn" polynomials. The 
recurrence coefficients can be computed safely from a suitable algorithmic use of 
the Shohat- Freud equations [LeQ], or from asymptotic expansions [CIS] (see [Nev2] 
pp.462-463). 

4. The same orthogonal polynomials reappear in special solutions of important dif- 
ferential equations of mathematical physics ([Bes] [KvMj; Shohat- Freud's equations 
are called "discrete Painleve equations" in [Fok2]), as well as in continued fraction 
expansions of special irrational numbers [Chu2]. 
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Any advance in one of these fields is liable to benefit to the other ones, although the 
dialog is not always obvious: established theorems may sometimes have poor constructive 
contents and be unable to inspire valuable algorithms; explicit formulas (using for instance 
exotic special functions or high-order determinants) may be delightful solutions for some 
people and useless nightmares for other ones; successful numerical tricks or self-consistent 
"vlnsatee" may be out of reach of contemporary methods of proof, etc. 

Let us return now to the problem considered here: to deduce properties of the 
recurrence coeflBcients a„ and bn from the measure dfj,{x). The Chebyshev orthogonal 
polynomials are related to measures involving the square root of a polynomial of degree 
2 and yield constant recurrence coefficients (the simplest case!). The classical orthogonal 
polynomials (Jacobi polynomials and their limit cases) have a known measure and known 
recurrence coefficients (a^ and are special rational functions of n). We may define a more 
general class by allowing and bn to be general rational functions of n (PoUaczek class 
[Chi]) but then the orthogonality measure becomes difficult to control. Natural extensions 
of Chebyshev polynomials are related to measures involving the square root of a polynomial 
of degree > 2. One finds then an oscillatory behaviour of the recurrence coefficients ([GVl] 
[GV2] [Gr] [I] [VA]), whose description may even need elliptic or hyperelliptic functions 
[Ak] [Apt] [GaN]. We will deal here with a further extension, the so-called semi-classical 
class (to be defined in the next section) which represents still a reasonable trade-ofi' between 
measure description (easy and containing useful cases) and the possibility of description 
of recurrence coefficients (interesting nonlinear relations). Much of the work was already 
done in the end of the nineteenth century by Laguerre [Lag] who recognized (in 1885) that 
special cases (worked by Jacobi and Borchardt) would indeed involve elliptic functions. 
As he could not establish the general recurrence coefficients behaviour, we could suspect 
that special functions still unknown in the nineteenth century would be needed. . . Painleve 
transcendents will indeed appear, and they were investigated in the early twentieth century 
(see the foreword of [Pain]). 



2. Formal semi-classical orthogonal polynomials. 

Orthogonal polynomials pn are usually defined through a measure, so to satisfy 
(1). The construction of these polynomials only needs the sequence of moments //fc = 
Jg x''dn{x), k = 0, 1, . . . Formal orthogonal polynomials are only related to a numerical 
(real or complex) sequence of numbers fik^k = 0, 1, . . ., ignoring whether these numbers 
are actual moments of some weight or distribution on some support or not. The 
polynomial Pn{z) — ■jnz'^ + ln,iz^~^ + • • • + Pn(0) is then obtained from the equations 

InfJ'n+k + 7n,lMn+fe-l H \- Pn{0)lJ,k = 0, A; = 0, 1, . . . , 77, - 1 and 7n[7nA*2n + 7n,lA*2n-l + 

■ ■ ■ + Pn(0)/Un] = 1- These equations can be solved for any n = 0, 1, . . . if the Hankel 
determinants built with ^uq, . . . , H2n do not vanish ([Bre], [dBvR] § 7 , see the definition of 
regular formal orthogonal polynomials on p. 47 of [Dra] § 1.1-1.3). 

If we define a linear form C on the space of polynomials by C{x^) = n = 0, 1, . . ., 
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the polynomials Pn satisfy jC(pnPm) = ^n,mj'iT',m = 0,1,... ([Mar], [Mar2], where C is 
written jCq). 



Regular formal orthogonal polynomials always satisfy the recurrence relation (2), 
with po = 70 = aiPiiz) = {z - 6o)po(^), and where ho = -71,1/71, 

a-n = 7n-i/7n, K = 7n,i/7n " 7n+i,i/7n+i, n = 1, 2, . . . ([Bre] , [Dra] § 1.4). 



By introducing the formal series 

oo 

f{z) = J2l^kZ-''-\ 



(3) 



the equations for pn are summarized as 

f{z)Pn{z) ^pi%{z)+er.{z), s^iz) ^ + 0{z—^) , (4) 

where p^^-i ^ polynomial of degree n — 1, (associated polynomial to Pn). These 
polynomials, as well as the e^'s, satisfy the same recurrence relations (2), but withp^^ = 0, 
Pq'^ = //o7i = 1/(qi7o)- The following relation 



PnPn-2 - Pn-lPn-l = Pn-l^n - Pn^n-1 = -l/^r 



(5) 



is well known ([Chi] , [Frl] , etc.) From the recurrence relations (2), we have the main 
terms in the expansions of p„ and £„, which will be useful later: 



Pn{z) = 7^ 



£n{z)^-fn^ 



n-1 



J2 bibj-J2^ 



2 1 ^ 



.i<j<n 



-n-2 



n+1 

i<j<n 



2 I ^ 



(6) 



(7) 



(for the latter one, use -fn£n{z) = {z - bn) ^-fn-i£n-i{z) + {z - bn) ^a^.^-fn+iSn+iiz)). 



Of course, if we happen to know a true function of the complex variable z having the 
asymptotic expansion (3) when ^ — > oo in some way, and if this function is analytic outside 
a set S made of contours and arcs, we may use a Cauchy-like integral representation 

f{z) = [ w{x){z-x)-^dx, z^S (8) 
Js 

allowing to recover the convenient description in terms of a "weight function" w, but the 
description is not unique and w may be complex. We then have an integral representation 
of the form C: C(p — Jg (p{x)w{x)dx. For instance the Bessel orthogonal polynomials are 
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defined by /x„ = l/n!,n = 0,1,... and can be considered as orthogonal with respect to 

the complex weight (27rz)~^ cxp on any contour containing the origin in its interior. 
Remark that orthogonality of two complex function ip and always involves here the 
product and not the product Lpi\) (as in [StT]). 

For an example showing how formal orthogonal polynomials can be investigated 
through their generating function of formal moments (3), consider f{z) — \A{z) — 
B{zYl'^\IC{z), where A,B and C are given polynomials {formes du second degre in 
[Mar2] p. 122, Def. 7.4). Such a function can be represented as (8) outside a systems 
of cuts S joining the zeros of B in some way. Here w{x) will have the form w{x) = 
±(7ri)~^i?(a;)^/^/C(a;) ([N]§1.2 & 4.3.1). If B has only real zeros, this is a way to introduce 
special orthogonal polynomials on several intervals (the intervals where B{x) < 0). Now, 
(4) gives here —B^/'^p^ = Qn + Csn, with g„ = —Ap^ + Cpl^}_i. Squaring yields 
~ — where must be a polynomial of degree bounded by a constant, as 
the left-hand side is a polynomial, and as the right-hand side is 2qnCen + C^^n- So, pn 
is such that the square of this polynomial times a given polynomial B equals the square 
of another polynomial plus a polynomial of bounded degree. This is enough for experts 
to describe Pn in terms of (hyper) elliptic function and integrals, theta functions, etc. (see 
[Ak] §53, [Apt], [Brez] pp. 296-298, [N]§4.3 , [Peh] ), and to discuss periodic features in the 
sequence of the recurrence coefficients ([GVl], [GV2], [Gr], [I], [Pehl] ). For arithmetic 
continued fractions connected to Pell's equation, see [Brez] pp. 39.43. 

A similar technique will now be applied to a more general class of functions /. 

Many special families of orthogonal polynomials have been studied. In most cases, 
the knowledge of a special family is considered satisfactory when an explicit formula for 
the recurrence coefficients and bn as functions of n is associated to a definite formula 
for the weight w, or measure of orthogonality, see for instance the final tables of Chihara's 
book [Chi] , whereas the starting point of the study may be generating functions, Rodrigues 
formulas, special functions identities, differential equations, etc. 

The simplest way to start the study of the class of semi- classical orthogonal 
polynomials is to define them through a differential equation of their function /: 

Definition: The sequence {pn{z) = 'JnZ^ + ■ ■ "j^^Lo is a set of formal semi-classical 
orthogonal polynomials if (3) holds with a function / satisfying the first order linear 
differential equation 

Wf = 2Vf + U (9) 
where W, V and U are polynomials {W ^ 0). 

This is equivalent to the existence of a linear recurrence relation of the form 
X^^=o("'^fe Vk)l^n+k — for the formal moments Hn [BeR]. 
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Moreover, only regular semi-classical orthogonal polynomials will be considered here, 
so that 7„ 0, n = 0, 1, . . . 

Of course, (4) must be possible with an expansion of the form (3), so that degree 
U < max( degree W — 2, degree V — 1). AH the classical families are recovered when 
degrees of W and V < 2 and 1. 

We will consider especially 

Definition: Generic semi-classical orthogonal polynomials are semi-classical orthog- 
onal polynomials where m = degree W > 2, degree V < m, the zeros Xi, X2, . . .Xm of W 
are distinct, and the residues ak — 2V{xk)/W' {xk) are not integers, k — 1,2, . . . ,m. The 
Jacobi polynomials correspond to m = 2. 

We have then: 

Proposition: Generic semi-classical orthogonal polynomials are orthogonal with 
respect to a (possibly complex) generalized Jacobi weight function 
w{x) = Ajir[^{x — Xk)°'^ on arcs Sj, j = 1,2, .. .m of the complex plane. 

Indeed, (9) has exactly one holomorphic solution fj{z) = cj^ + Cj^i{z — Xj) + ■ • ■ in a, 
neighbourhood of the singular point Xj, as the equations for the Cj/s are 2V{xj)cj^o + 
U{xj) = and W {xj)icj^i + ■■■ = 2V{xj)cj^i + ■■■,i — 1,2,... have exactly one 
solution, as V{xj) ^ and W'{xj)i — 2V{xj) = W'{xj){i — aj) cannot vanish (this 
can also be seen as a most elementary application of L.Puchs theory of linear differential 
equations). As 11 (a; — x^)'^'' is a solution of the homogeneous equation (9), one has 
f{z) = fj{z) + BjIlY^{z — Xfc)"'" near Xj, on one side of the cut. A Cauchy integral 
expression of f{z) will, after a distorsion of the integration contour (as in [N] §1.2), involve 
the difference of the limit functions /_|_ and /_ which is a multiple of 11^ (a; — a; /-)"*= on a 
cut. This gives w on S. Let w{z) be a continuation of w on some side of the cut, then we 
have 

f{z) = fj{z) + CM^) (10) 

near Xj. 

Non generic semi-classical orthogonal polynomials can be considered as limit cases, 
for instance, a weight expP(x), where P is a polynomial, is the limit of (1 -|- P{x)/N)^ 
when A?" — > oo ... See [Al] and [Bel] for other proofs and examples. 

Anyhow, as f+ and /_ along the two sides of a system of cuts are solutions of the 
same equation (9) , their difference must be a solution of the homogeneous equation: semi- 
classical orthogonal polynomials are orthogonal with respect to a {possibly complex) weight 
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function w satisfying 



Ww' = 2Vw 



(11) 



on a system of cuts, masspoints may also be present if / has poles. Examples have been 
given in [BoN], [HvRl], [HvR2] and [Sho]; the whole class of true positive semi-classical 
measures on real sets is given in [BLN] . 



Conversely, Shohat [Sho] develops the theory starting from a weight function 
satisfying (11) on an interval. Let us generalize this to a given set of arcs S, and show 
that (9) is recovered: if needed, we multiply W and V by common factors in order to 
have \i'mW{x)w{x) = when x tends to any endpoint (eq. (6) of [Sho]). Then, from (8), 
W{z)f{z) = l^(a;)t(;(a;) (2;— a;)~^cia; plus a polynomial {Jg[{W{z)—W{x))/{z—x)]w{x)dx 
is a polynomial in z). The derivative gives 

{W{z)f{z)y = - I W{x)w{x){z-x)-^dx+ pol. = I {W{x)w{x))'{z-x)-^dx+ pol. , 
Js JS 

by integration by parts, using Ww ^ at the endpoints of S. As {Ww)' = iW' + 2V)w, 
and fg{W'{x) + 2V{x)){z - x)~^dx = {W'{z) + 2V{z)) Jg{z - x)~^dx+ a polynomial, we 
find indeed Wf — 2Vf+ a polynomial, i.e, (9). 



3. Differential relations and equations for formal semi-classical orthogonal 
polynomials. 



Now, we go further, following Laguerre ([Lag] sec. 2, see also [HvRl], [Per] § 76): 
from (4) and (9), 



O^W 



Pn-l _|_ fn 
Vn Pn 



2V 



W\p^^\'pn-p'nP'nU] ' ^VPn^lPn ' Up. 

pI 



Pn-l _|_ 
Pn Pn 

(1) 



u 



+ W 



Pn 



-2V- 



Pr. 



SO, 

G„ = W-bi-l -PnP^n-l\ ' '^Vp^n-lPn ' Up. 

is a polynomial of degree bounded by a constant, as 



(12) 



©n = -pIW 



Pr 



+ 2VenPn = W{enp'n - s'nPn) + ^VSnPn 



(13) 



is bounded by a power < max( degree W —2, degree V —1) for large argument. For given 
W and V, (13) with (6) and (7) allow to give Qn in terms of n and the recurrence coefficients 
a's and 6's. Moreover, expanding (13) up to negative powers of z yields equations for these 
coefficients. This is a first hint towards identities (Laguerre-Preud's equations) for the 
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recurrence coefficients of semi-classical orthogonal polynomials. See [BeR] for a technique 
involving Turan determinants. 



Identities like (12) involving orthogonal polynomials of arbitrary high degree on 
one side and polynomials of bounded degree with respect to n on the other side occur 
whenever one has a functional equation P{f) = for /, provided the elimination of 

/ in P{p:^ljpn + 0(^-2"-^)) = is simple enough. This happens if P applied to a 
rational function ip/il^ produces another rational function with denominator ^ of degree 
not much larger that twice the degree of '4>. Then, multiplication by this denominator 
^ will produce polynomials and, roughly speaking, products of ^ and the error term 
Q^^-2n-i-^ which will keep a small rate of growth at oo. Exemples of valid functionals 
P are quadratic polynomials (discussed in the preceding section: f — {A — B^/'^)/C =^ 
(C/ — AY — B = 0), linear differential operators of first order discussed here, both 
giving ^ = and Riccati differential operators (theory of Laguerre-Hahn orthogonal 
polynomials [Magi]). Difference operators may also be considered, they can leave things 
like ^{z) = pn{z)pn{z + h),^{z) = Pn{z)pn{qz), ctc. [Mag4] 



In the generic case, let W (z) = Il^{z — Xk) = z"^ — ^fc)-^"^ ^ + ■ ■ then 2V{z) = 

Wiz) j:^{a,/{z-Xk)) = (E^ «fc)^"^-' + [ErK^fc)-(Er^fe)(Er <^k)]z^-^ + - ■ •, using 
(13), (6) and (7): 

en{z)= (2n + l + ^«fcU'^-2 + 



+ 



n-l 



2n + 1 + ^ afc h>n - ^ + 2 ^ 6i + 6n + y^X^kXk) 



(14) 



From (5), replace 6^ in (12) by (p„_ipi^2i - Pnpi^22)«n©n: 

PnlliWpn' + Vpn + ttnOnPn-l] = Pn[Wp^^\' - Vp^^l^ + anOnPn-2 " ^Pn], which mUSt 

therefore have the form QnPnPn-u where Qn is a new auxiliary polynomial of bounded 
degree. Using again (5), one has 

fin = Onl^bi-lVn-l - PnPn-2] " «nl^bn-lPn-l + PnPn-2\ " anUp^p^^i 

And, with (4): 

fin = anW{en-lPn - ^nPn-l) + anV{en-lPn + £nPn-l)- (15) 

This yields the two differential relations: 

Wp'n = {fin - V)pn - ttnOnPn-l 
Wp^^l,' = i^ln + V)p^^l, - anQnp'nU + Up^ 
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We get rid of the Upn term of the second equation by forming an equation for using (9), 
and subtracting the second equation: We'^ = {Qn + V)£n — ttn^n£n-i- We recover the form 
of the first equation by using (11): W{en/wy = (0„ — V)£n/'w — anOn^n-i/w. In order 
to have a differential system, we have to give y'n^i {y = p or e/w) in terms of yn and y-n-i- 
As yn satisfies the recurrence relations (2), Wy'^_i = (fin-i — V)yn-i — an-iQn-iyn-2 
turns easily as Wy'„_]^ = anOn-iyn + {^n-i - V -{z- bn-i)Qn-i)yn-i- Prom (15), (13) 
and (2), 

nn+i{z) = {z-bn)<dniz) -flniz), (16) 

SO we finally have the differential system: 



Y' = AY : 



Pn £n/u! 
Pn-1 £n-l/w 



1 

W 



Pn En/w 
Pn-1 £n-l/w 



(17) 



This differential system gives the whole differential history of the semi-classical 
orthogonal polynomials. Laguerre [Lag] and many other people ([AtE] [Hal] [Ha2] [Nev] 
[Sho] etc. ) have preferred the scalar second order form obtained from eliminating y^-i in 
Wy'n = i^n - V)yn - an^nyn-i and Wy'^_^ = anQn-iyn - iP'n + V)yn-i: 



(18) 



with Kn = (On - V)'Qn ' {^n ' V)&^ + Oni^l - V - 6, Gn-i )/ W^, which is a 
polynomial, as putting an+iy'n+i = - hn)y'n + yn - any'n-i (derivative of (2)) in 
On+iW^y^i+i = an+i{fln+i - V)yn+i - al^-^Qn+Wn, using again (2), and the differential 
equation (17) for Wy'^_i gives an expression of the form Ayn = Byn-i, with B = from 
(16), whence A = 0, which is 



{Z - bn){^n+l - ^n) = + On+l^n+l " On®n-l- 

Multiplying by (16) and summing on n, one finds 



(19) 



n-l 



(20) 



knowing that Qq = V. 



With Zn = {Ww/QnV^'^ym wc havc a form without first derivative 



le:: iQ'^w' + 2VLn - w" w" + 2n; Eo' 



MVen/ 2G„ 2Gn W AW^ 
used by R. Fuchs [RFu] in the case m = degree = 3. 
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2W 



W 



(21) 



Laguerre ([Lag], see also [GaN]) finds equations for the recurrence coefficients and 
the coefficients of and by using (16) and (19), keeping the degrees of and 
Qji bounded when n increases. We may express everything in terms of the recurrence 
coefficients alone, then the expansion of Qn, constructed on the same lines as (14), will be 
useful: 



nn{z)= [n+(J]"afe)/2] 

n-l 

— \'m /\ — — \rn ^ — \ 

~ ^Z^i + [1^1 (o^kXk) - (2^^ a;fc)(2^^ Qjfc) j /2 



n— 1 n— 1 /n— 1 \ 



+ 



,k<e<m 



. . (22) 



Consider for instance the case m = 3 (simplest generalized Jacobi polynomials): 
from (14), Qn is a polynomial of degree 1 with a known coefficient of z and a constant 
coefficient depending on the 6j's up to bn', from (22), is a polynomial of degree 2 with 
a known coefficient of z'^ and two other coefficients depending on the 6i's and the a^'s up 
to the index n — 1 (see example 1 in section 5) . The constant coefficients of (16) and (20) 
give nonlinear relations for a„ and bn- The meaning of the solutions of these recurrence 
relations for the recurrence coefficients of (2) is not obvious. Even the simplest relations 
found in nongeneric cases (as (30) or (31)) are baffling. 



The explanation in terms of Painleve transcendents and similar functions, i.e., 
solutions of remarkable high-order nonlinear differential equations in terms of a well-chosen 
parameter, will be given now. The derivation is based on the isomonodromy properties of 
(18). Later on, examples will show that a more elementary derivation is possible. 



4. Monodromy matrices and isomonodromy identities. 



D. & G. Chudnovsky remarked ([Chu2], see also (5.1.18) in [N]) how (18) has a 
form already investigated in the period 1890-1910 by authors working on isomonodromy 
deformations ([RFu], [Pain]). 



Let Y[z) be a fundamental matrix of solutions of the differential system Y'[z) = 
A{z)Y{z), defined outside a system of cuts joining the singular points (poles of A) of the 
equation. When z follows a contour about a singular point xj, let us solve Z'{z) = A{z)Z{z) 
with the initial value Z{zo) = ^^(^o) at a starting point on the contour. As long as no cut 
is crossed, Z{z) = Y{z). This is no more true when one or several cuts are crossed but. 
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when we come back in a neighbourhood of zq, the columns of the matrix of solutions Z{z) 
must be fixed combinations of the columns of the initial fundamental matrix of solutions: 
Z{z) = Y{z)Mj. This matrix Mj is called the monodromy matrix of Y' = AY at the 
singular point Xj (only regular singularities are considered here). 



Theorem 1. Generic formal semi-classical orthogonal polynomials satisfy differen- 
tial systems (1 7) with monodromy matrices 

M- ^•'■[■^ ~ exp(-27riaj)] 

■' exp{—2niaj) 

at the singular points Xj, j = 1,2, . . .m. 



Indeed, Pn and Pn-i are not modified after a circle about xj, but £„, Sn-i and 
w have a branchpoint there. According to the discussion made in the proof of (10), 
f{z) = fj{z) + Bjll'^{z — Xk)'^'' with some determination of the powers near Xj, near 
a side of a cut. By following a contour about Xj, fj returns to its previous value, 
but IlY^{z — Xk)"'' has been multiplied by exp(27riaj). The same happens with w. 

Therefore, from (10), Sn/w — {fpn — Pn-i)/w — {fjPn — Pn-i)/^ + CjPn becomes 

ex.^{-2'niaj){fjPn -Pn-i)/w + CjPn = exp(-27riQ;j)£„/w + [1 - ex.^{-2'niaj)]CjPn- 

This shows that the monodromy matrices of (17) at the singular points remain 
unchanged if the exponents remain unchanged and if the weight w on S is adapted 
so that the multipliers Ck remain the same. However, one may vary the positions of 
the singular points Xk- The quantities /, p„, an, bn, ©n etc. will then be subject to 
extrememely interesting isomonodromy deformations. Here is a sketch ([LD] HI, from 
p. 128 onwards), applied to the specific equation (17): 



Let the Xk depend on a single parameter t, and let us define the matrix 

dt ' 

as dMj/dt = 0, H does not change when z achieves a contour about Xj. So, H has no 
branchpoints at the x/s. To get a better view of what happens at the singular points, we 
expand H (using dety = l/{anw), from (5)): 



TT _ ^ Pn^n-1 - Pn-l^n + Pn-l^nW/w -Pn^n + Pn^n 

[Pn-ien-l -pn-lSn-l + Pn-l£n-lW/w -pn-lEn + Pn£n-1 - Pn£n-lW/w j 

(23) 

where the dot derivative is d/dt. From (4) and (10), one has = Snj + Cjwpn near Xj, 
where e^j is regular near Xj. The singular terms cancel nicely in the combinations of (23) 
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(remember that Cj = 0!); the ratio w/w has a simple pole at xj with residue —OijXj (as 
dj = 0). We are left with 

m 

H = H^ + J2Hj{z-Xj)-\ 

3 = 1 



with 



Pn—l^n,j Pn^n.j 
Pn—l^n—l,j Pn^n — l,j 



j = l,...,m, 



where the PrSgjn are the values at z = Xj. As W{xj) = 0, (13) tells that 0„ = 2Ven,jPn 
at z — Xj, and (15) with (5) gives — V + 2anVenjPn-i = ~V + '^O'n.V^n-ijPn Sit Xj. 
With aj = 2V/W' at xj, one finds from (17): 



A direct inspection of (23) when 2; — > oo gives, using (6) and (7), 
in the generic case, as w/w = — OikXk/{z — Xk) when z — > 00. 



inhn 

-in-l/ln-l 



Finally, the differential equations in t appear by working 

d'^Y/dzdt = dY/dz = {HY)' = H'Y + HY' = H'Y + HAY = 
= d^Y/dtdz = dY'/dt = (AY) ^ AY + AY ^ AY + AHY, 

whence 

A = H' + HA- AH. (24) 

C'etaient les cieux ouverts 
Stendhal 

This equation (24) has an incredibly inspiring form, explaining how this theory is 
related to integrable Hamiltonians, Biicklund transformations, Lax pairs, Toda lattices, 
solitons, etc. whereas the connection with orthogonal polynomials, special functions, 
continued fractions, Diophantine approximations has been worked with great virtuosity 
by G. & D. Chudnovsky [Chua] [Chub] [ChuO] [Chul] [Chu2],. . . 

In the generic case, we have for the residue matrices 

k=m . 

Aj — HqqAj — AjHqq + ^ ^ — [Af^Aj — AjAf-), j = 1, . . . , m 

fc=i ~ 
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called the Schlesinger equations (see [Chua]). 



Now, we show how these equations lead to differential equations for the recurrence 
coefficients: 



Theorem 2. Let an and bn be the recurrence coefRcients of (2) for generalized Jacobi 
orthogonal polynomials related to a ( possibly complex ) weight of the form IVP'(x — Xj)^^ , 
on a set of arcs joining the xj 's, where at least one of the Xj 's depend on a parameter t. 
Then, we have the Toda equations 



a- 



On _ 1 Y> (6>n(a^fc) - Qn-l{Xk))Xk 



(25) 



k=l 



W'{xk) 



(26) 



where W{x) = Il'^{x — x^), and 6^ and Vln are polynomials introduced in (12) — (15). 



Indeed, from (17), the residue matrix Aj is 



^niXj) - V{Xj) -anOn{Xj) 



W'(Xj) [ anQn-l{Xj) -ilniXj) -V{Xj) 



we have 



H A — A H 



7n/7n + 7n-l/7n-l 
W'{X,) 



@n{Xj) 
®n-l{Xj) 



X 



X 



W'(xj)W'ixk) 

an{OniXj)On-liXk) " Qn{Xk)Qn-l{Xj)) 2 (O,, (,Xfc )n„ (,Xj ) - OriiXj)^niXk)) 

2{Qri-l{xk)^n{Xj) - Qn-l{Xj)^n{xk)) an{Qn{Xk)Qn-l{Xj) - Qn{Xj)Qn-l{xk)) 



The Schlesinger equations for the off-diagonal elements of Aj are 



dt W'{xj) 



Xj Xk ©n ('^A; )^n ('^j ) ®n(-^j)^n('^fc) 



lnW'{Xj) t^,Xj-Xk 



dt W'{Xj) 



W'{xj)W'{xk) 

Xj - Xk Qn-l{xk)^n{Xj) " Bn-l{Xj)VLn{xk) 



(27) 



7n-i W'{xj) trrxj-xk 



W'{xj)W'{xk) 



where a^7„ = 7n-i =^ Cin/O'n + in/in — in-i/ln-1 has been used. Increasing n by 1 in 
the second equation and adding to the first one. 



= -4 



7ti {,X 



In W'ix 



. Qn{Xj)Qn{Xk) 
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W'{xj)W'{xk) 



where (16) has been used. At this point, we don't have to avoid the term k = j 
anymore in the sum. Moreover, as any polynomial P{z) = tiqz'^~^ + ■ ■ ■ satisfies ttq = 
YJ^P{xk)/W\xk), (coefficient of ^"^ in P{z)/W{z) = P{xk) / {{W {xu){z - xu))), 
and as the degree of is m — 2 ((14)), Xj disappears from the sum: 



ln~ 2^ W'{xk) 

and (25) follows from a„7„ = 7n-i- 

Now, we come to the first diagonal element of the Schlesinger's equations: 

d Q.n{Xj) - V{.Xj) ^ 2 \p - ^n{Xj)Qn-l{.Xk) - Qn{xk)Qn-l{Xj) , . 

dt w'{x,) w'{x,)w'{xk) ' ^ ^ 

f . , Qn{x)Qn-l{y) - ^n{y)Qn-l{x) . i -l 

tor J — l,...,m. As is some polynomial, say 

x-y 

gTp^gX^y'^, of degree max(p, q) < m — 2 in x and y, the term k = j may be included 
in the sum as before. Still using Pi^k) /W' {xk) — for polynomials P of degree less 
than m — 1, Xj may also be removed, and 

^ 1 d Q.n{Xj) -V{Xj) ^ 2 Xk q 

^ _ ^ . at W'ixj) - l^'(a;,) ^ ^^'^^^ ^ (^ - a;,) W^'(a;,) 

m . „ 
2 V ^ Xk 



_ 2 Y^ -^fc Q 
m 

_ 2 ^ 



i^fc Qn(^)Qn-l(a:fc) - On{Xk)Qn-l{z) 

~[W'{xk) W{z){z-Xk) 



The left-hand side is 

d / 1 Cln{xj) — V{xj) \ 9 / 1 \ f2„(a;j) — F(a;j) 



dt \fr{z-Xj W'{xj) j fr[dt \z-XjJ W'{xj 

_ d_ ilnjz) - V{z) _ ^ Xk ^njXk) - VjXk) 

~ dt W{z) ^^{z-Xkf W'{xk 

and we take the z~'^ term in the expansion about oo, the right-hand side vanishes as degree 
0„ < m — 1, and using (22) in the left-hand side: 

yields (26). 
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In concrete situations, (25) and (26) will be used, together with other non differential 
identities (Freud Laguerre equations for the recurrence coefficients), but we may prefer to 
return to (24), or even use ad hoc differential relations. In the generic case, (27) and 
(28) for J = 1, 2, . . . , m give a differential system for 2m unknowns Qn{xj) and Qn{xj), 
j = 1,2, ... ,m, when the quantities a'^Qn-iixj) are eliminated with the help of (20) at 
Xj (recall that W{xj) = 0). However, considering from (14) and (22) that there are only 
2m — 3 unknown coefficients in Gn and O^, further eliminations are possible. We start 
with an example of generic semi-classical weight with m = 3. 



5. Example 1. Generalized Jacobi weight with three factors 

So, W{z) = z{z -l){z- t), V{z) = {az{z - t) + (3{z - l){z - t) + -fz{z - l))/2, the 
support S joins 0, 1, and t in some way, or is an arc joining only two of these points. (14) 
and (22) yield readily 



with Vn = 2n + 1 + a + /? + 7, 'dn = Vn{K - I - t) + 2^^^ ^ hi -\- hn + a + ^t, 

«n = Eo"' - ('^n - 1)(1 + t)/2 + (a + 7t)/2, and a;„ = Eo"'(^' - + 1)^* + 2a?) - 
{t + l)(a + 7t)/2 + - l)t/2 + (a + 7^^)72 + v^a^ 



(25) and (26) are here, with xu = Sk,3, W'ixs) = W'{t) = t{t - 1), 

On ^ -2 + {Un + l)bn - (^^n " 3)6n-l ; ^ bnjbn - 1) + jl^n + 2)a^+i - {u^ - 2)a^ 

an ~ 2t{t - 1) ' ~ ^(^ - 1) 

One would have a true differential system if bn-i and a^^j^ were simple functions of bn and 
ttn, but this does not seem to be the case here. So, we try with the unknowns k„ and 
Un instead. In (27), using [e„(x)0„(y) - e„(y)0„(x)]/(y - x) = (z/„ - l)z/na:^y/2 + ^?n('^n - 
l){x + y)/2 + Cn, withCn = ni^n ~ ^n^ni with Xj — 0, 1, t. One finds three equations which 
are all equivalent to 

t{t-l) • 



In (28), using [Qn{x)On-i{y) - Qn{y)Qn-i{x)]/ {x - y) = Un'dn-i-'^nt^n-i, One finds two 
independent equations 

(^n al{l^n'^n-l - '^nl^n-l) 



UJr. 



t t{t-l) 



«n = XT: 7T + -. T + 



2{t-l) t-1 t{t-l)' 
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Now, the three non differential equations (20) at a; = 0,1, t allow to eliminate a^dn-i-. 



, u;l-pHy4 _ (K-l)/2 + /.. + a;„)^-a^(t-l)V4 

an^n-l — a ' — — Cb^l^n-l, 

and a third equation allowing to eliminate either nn or uj^i actually it is simpler to 
give everything in function of Qn'- from a^_][dn-i/t - {'dn-i + i'n-i)/{t - 1) + {'&n-i + 
Un-lt)/{t{t - 1))] = 0, 

a^Mt - l)/4 ^ PH/A ^^^nt{t - l)/4 



K - l)7^n[z^nt(t + 1) + + t+l)]/4+ [Vnt + ^n(t + l)]Cn + Cl/ {^n ' 1) 

This allows to give as a function of and C^, so to complete (29): 
~ H ^t: — 1- r + — 7T : 1" 



t{t-l) 2(t-l) t-1 t{t-l) t ^(^-1) 

Using the preceding calculations, a^'&n-i and a^fn-i are replaced in terms of Hu and Um 
then Kn = (Cn + Vn'jJn)/'&n is used, and a;^ is finally replaced as a function of Cn; and what 
comes out is 



t{t-l)\ 4(iy„+^n) ^^n 



whence, at last, with Cn — '*^n('*^n + l)/2 = t(t — l)'i?n/2: 



2 \iyn+^n l^nt + "^n J \t t-1 V^t + {}n 

aH^i^nt + l3n) /S^il^n + ^n){^ut + K) ^ ( 1 - )^n (^^n + ^n) ^ 



2t2 (t - 1) (zy^ + ^?„) 2t(t - 1)2^?^ 2t(t - l){Vr^t + 7?„) 



2t2(t-l)2 

which is a Painleve equation of the sixth kind ([In] § 14.4) in —"dn/i'n (the zero of 0^) 
([Chua] p. 399-402, explaining works of R.Fuchs on equations of form (21)). 
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We can return to a„ and 6„ as functions of and C,n by using again a^Vn-i as a 
function of 1?^, and uJn (and I'n-i =2?i + q; + /? + 7 — 1 is known) and taking from 
2Kn—'&n = — (22^^ — 1) (1+^) " ('^n + l)&n- Inverting the connection should give a (probably 
algebraic) differential system involving only and 6„ (will somebody do that?) 



6. Example 2. exp(x'^/3 + tx) on (x : x"^ < 0}. 



Much simpler identities occur when the weight w is the exponential of a polynomial, 
so that w' /w is a polynomial itself. Recall (end of Section 2) that W{x)w{x) — > when x 
tends to the endpoints (if any) of the support S. We want the simplest case {W{x) = 1), so 
that the support cannot have finite endpoints, but must end on directions where w{x) — > 0, 
with at least one complex direction (or else all the moments vanish). So, we can take the 
set {x : < 0}, or for instance only {x : aigx = ±27r/3}, or also some equivalent contour, 
as {x : Re a; = a positive constant} leading to Airy functions and integrals ([Chu2] [Marl]). 
The weight can be considered as a confluent generalized Jacobi weight with singular points 
at oo: w{z) = limAr^cx)[l + (^^/3 + tz)/N]^, with an exponent N independent of the 
parameter t. As (25) and (26) hold for distinct finite singular points, we return to (24) 
assumed to be still valid: here, W{z) = 1, '2V{z) = w'{z)/w{z) = z'^ + t. Working (13) 
and (15) about oo, we have 



Quiz) = Z + br. 



^n{z) = {z^ + t)/2 + a: 



2 



Pushing (13) and (15) up to the z~ 
equations, i.e., the identities 



term, one finds the corresponding Laguerre-Freud 



a. 



+ «n+l + ^n+* = 0, n+al{bn + bn-l)=0. 



(30) 



We compute H in (23) up to the 0(1) term, as H is now expected to be a polynomial (see 
[Fed] § 2), taking care oiw/w = z: 



A = 



a: 



an{z + bn-l) 



i{z + bn) 



-az 



H = 



On 



-in-l/ln- 



The diagonal elements of (24) yield 2d„ = a„(6„ — bn-i)-, and the off-diagonal elements: 

'^inhn+bn = 0, dn6„+a„6^ = an(7„/7n+7n_i/7n-i)6n-2a^-ani and d^^n-l+an^n-l = 
—aniin/ln + 7n-i /7n-i )^n-i + 2a^ + o^t. Uslug (30), all these equations are compatible 
with the differential system 




bn + 



h — —h^ — 2n^ — f 



n 

''n 



which is the differential system equivalent to the second Painleve equation for (— 6„,4a^) 
([Chu2], [Ge] p. 339). The connection with Painleve transcendents can lead to advances in 
the solution of the problem posed by Maroni in [Marl]: when do we have ai, a2, • • • 7^ in 
(30)? The problem is now to localize the zeros of solutions of special Painleve equations. 
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7. Example 3. exp(-a;^/4 - tx^) on R. 



This is the simplest nontrivial Freud's weight, and the corresponding orthogonal 
polynomials have been much worked ([BoN] [Fr2] [LeQ] [Lub] [Mag2] [Mag3] [NeV] 
[GFOPCF] [NeV2] [Sho] ). As for example 2, we expand (13) and (15) with W{z) = 1 and 
2V{z) = -z^ - 2tz: 



Qniz) = -z' -2t-ai- a^+i, Qniz) = -zy2 - (a^ + t)z. 

A relation between the a^'s is found by expanding (20), equating the z'^ terms gives 

al{al_^ + al + a^+i) + 2tal = n, n = 1, 2, . . . (ao = 0) (31) 

a relation which seems to have been found by Shohat ([Sho] ), rediscovered by Freud 
[Fr2] and Bessis [Bes]. Remark that we have a degree of freedom on oi: this is because 
the weight can have the real axis and the pure imaginary axis in its support, with 
w{x) = Aexp(— a:'^/4 — tx'^) on the pure imaginary axis, and the preceding results hold 
for any A, so ai is some function (which can be computed from first moments) of A. 
However, if it is requested that all the a^'s are positive, the solution is unique and can be 
computed efficiently ([LeQ], see also [Nev2] p. 470). Now, (17) and (23) are computed: 



A = 



—aiz 



-aniz'' +2t + al_^+al) 



ttniz'^ + 2t + al + al^i] 



.2 _^ 0+ I ^2 I ^2 

z^ + {al + 2t)z 



H = 



In — Cin 
-anZ 



QnZ 

in-\hn-\ + Z^ + al 



The equations from (24) amount to be equivalent to 



which, with a„7n = 7„_i, gives 



7n 
In 



(32) 



On 

an 



(33) 



Actually, (32) (and (33)) can be recovered by quite elementary means: let {pn{x]t)} 
be the polynomials orthonormal with respect to an even measure of the form da{x;t) = 
exp{—tx'^)da {x; 0) on some support S, we have then for the monic orthogonal polynomials 

Pn/ln- 



d 1 d f f Prr{x]ty ^ 



exp(— ta; ) da{x; 0) 



dt ll dt Js V 7n 

2 [ Pn{x]t)\ / + 2\ 7 / ri\ ^n'^^n+1 



2 ^2 _L ^2 

^2 I fn\ , } I gxp(— ta;^) dcr(a;; 0) ^ 

7n / 7n 
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using x^pn = anan-iPn-2 + ia'l. + a'l+i)Pn + an+ia'n+2Pn+2 from (2) when bn = 0, and that 
the derivative in t of a monic polynomial must be of degree < n. Conversely, it has been 
shown that (33) implies that the the coefficients of the recurrence of orthogonal 

polynomials with respect to a measure of the form exp{—tx'^)da{x; 0) where da{x; 0) does 
not depend on t [KvM] [Mo] (see also [Pral], [Fra2]), [Y]. 

It is even probably possible to recover the information given by (24) for all the 
semi-classical orthogonal polynomials by more elementary means, but the connection with 
monodromy theory, interesting on its own right, has more advantages: for instance, it 
is known that the differential equations produced by (24) have the Painleve property 
(foreword of [Pain], see [Mai] for a modern proof), i.e., movable singular points can only 
be poles (see [Cha], [In] chap. 14). No wonder that the classical Painleve transcendants 
appear in these examples. 

Now, we get an equation for the single a„ using (31): let tt„ = a^, from = 

Un{Un-l - Un+l), 

Un = Un{Un-l " Un+l) + Un{Un-l " ^n+l) 

= Un{Un-l - Un+lf + Un[Un-l{Un-2 " Un) - Un+l{Un - Un+2)] 
= Un{Un-l - Un+lf + Un[n - 1 - 2tUn-l - Un-l{Un-l + 2tt„) + 
+ n+l - 2tUn+l - Un+l(Un+l + 2Un)] 
= Un[2n - 2{Un + t){Un-l + Un+l) - 2Un-lUn+l] 

= Un[2n - 2{Un + t){Un-l + Un+l) - {Un-1 + Mn+l)^/2 + {Un-1 - Un+lf /2] 
= Un[2n + 2{Un + t)^ - {Un-1 + 2Un + Un+l + 2t)^ /2] + (Un)^ / {2Un) 
= Un[2n + 2{Un + t)^ - {n/Un + Unfl2\ + {Unfl{2Un) 



Un 
2 



\Un J 



+ 



2Ur 



+ (3^2 + 2tun - n) [ul + 2tun + n) , 

^Un ^Un 



which is a special case of the 4*"^ Painleve equation 



j/=| + ^+4ty^ + 2(t^-«)y+^ (34) 
with a = -n/2 and /3 = -nV2 [Bu] [Fokl] [Fok2] [Ge] [Ok]. 

For an = ^/u^i one has a form without first derivative: 

4a3 an = (3a^ + 2tal - n) (a^ + 2tal + n) . (35) 



Let an = n ^^^cin and t = n ^/^t, then we have another form 

4a^^a„ = n2(3a^ + 2ta^ - l)(a^ + 2ial + !)• (36) 
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What can be the use of these equations? To explore these things, let us first look at 
the graph of some On's computed with the Lew & Quarles method [LeQ]: 



ai,...,aio (only ai,...,04 arc marked) tend to be close to the zeros of the right- 
hand side of (36) (thick line). In particular, an(t) ~ l/\/2t when t — > +00: a\ = 
IJ'^I \J^Q = f^^x'^w{x)dx/ f^^w{x)dx, where w{x) = exp(— — tx^). From [Erd] 

p. 119, /Uo = V Tv\^exp{t'^ /2)D_i/2{tV2) (parabolic cylinder function). When t +00, 
jjQ ~ A/yiTt ([Erd] p.l22), so = ^2/1^-0 = -/io/A*o ~ l/(2t). From (31), if ai, 02, . . . , a^-i 
are 0(t~^/^), ~ i/n/ (2t) when t — >• 00. When t — >• —00, /xq ~ constant t~^/^exp(t^) 
([Erd] p. 123), so ai ~ \/—2t. The figure suggests that On(t) ~ V— 2t when t — > —00 and n 
is odd, while an(t) ~ •>/— l/(2t) when n is even. 

One of the most interesting uses of expressions of recurrence coefficients where n is 
not bound to be an integer is to define general associated orthogonal polynomials, 

i.e., polynomials defined by a^^^j^xV^nXxiA = -^n^v)Vn^ {z) - an^vV^n-\{^)'> ^^d degree 
'P'n^ = n (as in [AW], [ILVW]). So, let us define ai, as some solution of (35) with n replaced 
by u: 

4ala^ = {3al + 2tal - iy){at + 2tal + v), (37) 



v al dp 



2a2 2 * ^ a 



1/2 

satisfies the 



where is a given complex number. Then, y = 

same equation (37), but with v replaced by z/±l (Schlcsinger transformation, [Fokl] §3.3). 
Indeed, derivating y2+t-|-(a^ — z//a^)/2 = ^a^j/ajj yields 2yy = ^{i>±l — 2aly'^ — y'^ — 2ty'^) 
and a new derivation establishes the property. So, the definition makes sense and (31) still 
holds with u. There are still two degrees of freedom in (37), but they are removed when 
suitable boundary conditions are fixed ([DeCl] [DeC2]). Here, we just have to impose 
ttj^ = 0(t~-^/2) when t — > -|-cxo ([Yos], quoting Malmquist; the point being that ai,{t) must 
have an asymptotic series when t — ^ +00 for fixed z/, the relation with the dual situation, 
i.e., t fixed and u — > +00 is striking, see Section 4 of [Wi]). In summary: 
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For any real or complex u, the associated Freud orthogonal polynomials (which 
are related to the weight exp(— — tx'^ on M when u = 0) have recurrence coefficient 
ajyj^i{t), 01^+2 (^)) • • where a^{t) is completely defined as the solution of 

Aala^ = (3aJ + 2tal - ^i){al + 2tal + 
which remains 0(t~^/^) when t — > +00. 

For the associated polynomials themselves, we can now construct Qv+n and O^u+ni 
therefore a differential equation (18) with index v -\- n. Let ip^j^n and V'j^+n be two 
independent solutions of this differential equation (in z). Following Hahn ([Hal] eq. 
(17)), p^n"* = [fu+ni^iy-i - i^u+nfiy-i) / {fui^u-i " il^ufu-i)- It cau then be shown that 
fi, — limn^oo p^^^^^ / pn^ satisfies a Riccati equation (Laguerre-Hahn class [Magi]). 

8. Example 4. {x — t)P exp(—x'^) on [t, 00). 

The corresponding orthogonal polynomials are called (when t = p = 0) the 
Maxwell polynomials in [BeR], where other references can be found (p = 1: speed 
polynomials in [CIS]). This case is closely related to the preceding one: put x = 
t + in pn{x)pm{x){x — tY ex.-p{—x'^) dx = 5m,n to find that P2niu) = 

2-(p+i)/2 exp(— i^/2)p„(t + w^/2) is the orthonormal polynomial of degree 2n with re- 
spect to the weight w{u) — I'up''"'"^ exp(— — tu^) on M. So, we have a„ = a2nS2n-i/2 
and bn = t + {a\^ + ai^+i)/2 ([Chi] , etc.). 

For the a^'s, we still have a„ = a„(d^_j^ — a^_|_j^)/2 as before, but a slightly 
different recurrence relation a^(a^_i + + a^_|_i + 2t) = n + (2p + l)odd(n), where 
odd(n) = (1 - (-l)'')/2 [Fr2] [Mag2]. Working this yields now (u^ = dl). 

tt^ 3 3 /on ,,l + 3(-l)'^\ (n+ (2p + l)odd(n))2 

i.e., the Painleve 4*^ equation (34) with a = -n/2 - (2p + 1)(1 + 3(-l)'")/4 and 
P=-(n+ {2p + l)odd(n))V2. 

Many almost-classical orthogonal polynomials (see [Chin], [CIS] and references in 
[BeR] and [Gau]) could still be worked, and the simplest of them will likely be related 
to other Painleve transcendents (perhaps not the first one. . . although [Fok2] finds first 
Painleve transcendents as solutions of a limit case of (31)) At least a new case is briefiy 
presented now: 

9. Example 5. Beyond Painleve: exp(— — tx^) on M. 

With {un-1 - Un+i) stiU holds, but the recurrence relation is 

somewhat more complicated than before [Fr2] [Mag2] [Mag3]: tt„(w^_2«n-i + 'U'n-i + 
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2Un-lUn + ul + 2UnUn+l + Itn-l^n+l + ^l+l + «n+l^^n+2 + 2t) = n, forn = 1, 2, . . . As a 

first step, one has a differential system for Un-i, ■ ■ ■ ,'Un+2 by eliminating Un-2 and u^+s 
from the recurrence relation: 

iln-l = Un-lUn-2 " ^n-l^n = 

= — - 2t- u1_i - 3Un-lUn -ul- 2UnUn+l - Un-lUn+1 - " '"n+lWn+2, 

iln = Un{Un-l - ^n+l), 
Un+1 = Un+l{Un - Mn+2), 
^n,+2 = Un+lUn+2 — Un+2Un+3 = 

= \-2t + tt„+2 + «n^^n+2 + SUn+lUn+2 + U^+l + 2UnUn+l + + «n-l«n, 

which can still be transformed. . . This case in considered in [Fok2] . 
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